[INFO   ] Starting gmx_MMPBSA v1.5.2
[DEBUG  ] WDIR          : /media/romi/My Passport/November 3rd 2023/6umt-rs1-310k/free-energy
[DEBUG  ] AMBERHOME     : /home/romi/anaconda3/envs/gmxMMPBSA
[DEBUG  ] PYTHON EXE    : /home/romi/anaconda3/envs/gmxMMPBSA/bin/python
[DEBUG  ] PYTHON VERSION: 3.9.12 (main, Apr  5 2022, 06:56:58) [GCC 7.5.0]
[DEBUG  ] MPI           : /home/romi/anaconda3/envs/gmxMMPBSA/bin/mpirun
[DEBUG  ] ParmEd        : 3.4.3
[DEBUG  ] OS PLATFORM   : Linux-5.15.0-83-generic-x86_64-with-glibc2.31
[DEBUG  ] OS SYSTEM     : Linux
[DEBUG  ] OS VERSION    : #92~20.04.1-Ubuntu SMP Mon Aug 21 14:00:49 UTC 2023
[DEBUG  ] OS PROCESSOR  : x86_64

[INFO   ] Command-line
  gmx_MMPBSA -O -i mmpbsa.in -cs md_0_10.tpr -ci index.ndx -cg 18 19 -ct system_reduced.xtc -cp topol.top -o FINAL_RESULTS_MMPBSA.dat -eo FINAL_RESULTS_MMPBSA.csv

[DEBUG  ] |Input file:
[DEBUG  ] |--------------------------------------------------------------
[DEBUG  ] |Input file generated by gmx_MMPBSA (v1.5.2)
[DEBUG  ] |Be careful with the variables you modify, some can have severe consequences on the results you obtain.
[DEBUG  ] |
[DEBUG  ] |# General namelist variables
[DEBUG  ] |&general
[DEBUG  ] |  sys_name             = ""                                             # System name
[DEBUG  ] |  startframe           = 10                                             # First frame to analyze
[DEBUG  ] |  endframe             = 4990                                       # Last frame to analyze
[DEBUG  ] |  interval             = 95                                        # Number of frames between adjacent frames analyzed
[DEBUG  ] |  forcefields          = "oldff/leaprc.ff99SB,leaprc.gaff"              # Define the force field to build the Amber topology
[DEBUG  ] |  ions_parameters      = 1                                              # Define ions parameters to build the Amber topology
[DEBUG  ] |  PBRadii              = 3                                              # Define PBRadii to build amber topology from GROMACS files
[DEBUG  ] |  temperature          = 310                                        # Temperature
[DEBUG  ] |  qh_entropy           = 0                                              # Do quasi-harmonic calculation
[DEBUG  ] |  interaction_entropy  = 0                                              # Do Interaction Entropy calculation
[DEBUG  ] |  ie_segment           = 25                                             # Trajectory segment to calculate interaction entropy
[DEBUG  ] |  c2_entropy           = 0                                              # Do C2 Entropy calculation
[DEBUG  ] |  assign_chainID       = 0                                              # Assign chains ID
[DEBUG  ] |  exp_ki               = 0.0                                            # Experimental Ki in nM
[DEBUG  ] |  full_traj            = 0                                              # Print a full traj. AND the thread trajectories
[DEBUG  ] |  gmx_path             = ""                                             # Force to use this path to get GROMACS executable
[DEBUG  ] |  keep_files           = 2                                              # How many files to keep after successful completion
[DEBUG  ] |  netcdf               = 0                                              # Use NetCDF intermediate trajectories
[DEBUG  ] |  solvated_trajectory  = 1                                              # Define if it is necessary to cleanup the trajectories
[DEBUG  ] |  verbose              = 1                                              # How many energy terms to print in the final output
[DEBUG  ] |/
[DEBUG  ] |
[DEBUG  ] |# (AMBER) Generalized-Born namelist variables
[DEBUG  ] |&gb
[DEBUG  ] |  igb                  = 5                                              # GB model to use
[DEBUG  ] |  intdiel              = 1.0                                            # Internal dielectric constant for sander
[DEBUG  ] |  extdiel              = 78.5                                           # External dielectric constant for sander
[DEBUG  ] |  saltcon              = 0.0                                            # Salt concentration (M)
[DEBUG  ] |  surften              = 0.0072                                         # Surface tension
[DEBUG  ] |  surfoff              = 0.0                                            # Surface tension offset
[DEBUG  ] |  molsurf              = 0                                              # Use Connelly surface ('molsurf' program)
[DEBUG  ] |  msoffset             = 0.0                                            # Offset for molsurf calculation
[DEBUG  ] |  probe                = 1.4                                            # Solvent probe radius for surface area calc
[DEBUG  ] |  ifqnt                = 0                                              # Use QM on part of the system
[DEBUG  ] |  qm_theory            = ""                                             # Semi-empirical QM theory to use
[DEBUG  ] |  qm_residues          = ""                                             # Residues to treat with QM
[DEBUG  ] |  qmcharge_com         = 0                                              # Charge of QM region in complex
[DEBUG  ] |  qmcharge_lig         = 0                                              # Charge of QM region in ligand
[DEBUG  ] |  qmcharge_rec         = 0                                              # Charge of QM region in receptor
[DEBUG  ] |  qmcut                = 9999.0                                         # Cutoff in the QM region
[DEBUG  ] |  scfconv              = 1e-08                                          # Convergence criteria for the SCF calculation, in kcal/mol
[DEBUG  ] |  peptide_corr         = 0                                              # Apply MM correction to peptide linkages
[DEBUG  ] |  writepdb             = 1                                              # Write a PDB file of the selected QM region
[DEBUG  ] |  verbosity            = 0                                              # Controls the verbosity of QM/MM related output
[DEBUG  ] |  alpb                 = 0                                              # Use Analytical Linearized Poisson-Boltzmann (ALPB)
[DEBUG  ] |  arad_method          = 1                                              # Selected method to estimate the effective electrostatic size
[DEBUG  ] |/
[DEBUG  ] |
[DEBUG  ] |# (AMBER) Possion-Boltzmann namelist variables
[DEBUG  ] |&pb
[DEBUG  ] |  ipb                  = 2                                              # Dielectric model for PB
[DEBUG  ] |  inp                  = 2                                              # Nonpolar solvation method
[DEBUG  ] |  sander_apbs          = 0                                              # Use sander.APBS?
[DEBUG  ] |  indi                 = 1.0                                            # Internal dielectric constant
[DEBUG  ] |  exdi                 = 80.0                                           # External dielectric constant
[DEBUG  ] |  emem                 = 4.0                                            # Membrane dielectric constant
[DEBUG  ] |  smoothopt            = 1                                              # Set up dielectric values for finite-difference grid edges that are located across the solute/solvent dielectric boundary
[DEBUG  ] |  istrng               = 0.0                                            # Ionic strength (M)
[DEBUG  ] |  radiopt              = 1                                              # Use optimized radii?
[DEBUG  ] |  prbrad               = 1.4                                            # Probe radius
[DEBUG  ] |  iprob                = 2.0                                            # Mobile ion probe radius (Angstroms) for ion accessible surface used to define the Stern layer
[DEBUG  ] |  sasopt               = 0                                              # Molecular surface in PB implict model
[DEBUG  ] |  arcres               = 0.25                                           # The resolution (Å) to compute solvent accessible arcs
[DEBUG  ] |  memopt               = 0                                              # Use PB optimization for membrane
[DEBUG  ] |  mprob                = 2.7                                            # Membrane probe radius in Å
[DEBUG  ] |  mthick               = 40.0                                           # Membrane thickness
[DEBUG  ] |  mctrdz               = 0.0                                            # Distance to offset membrane in Z direction
[DEBUG  ] |  poretype             = 1                                              # Use exclusion region for channel proteins
[DEBUG  ] |  npbopt               = 0                                              # Use NonLinear PB solver?
[DEBUG  ] |  solvopt              = 1                                              # Select iterative solver
[DEBUG  ] |  accept               = 0.001                                          # Sets the iteration convergence criterion (relative to the initial residue)
[DEBUG  ] |  linit                = 1000                                           # Number of SCF iterations
[DEBUG  ] |  fillratio            = 4.0                                            # Ratio between the longest dimension of the rectangular finite-difference grid and that of the solute
[DEBUG  ] |  scale                = 2.0                                            # 1/scale = grid spacing for the finite difference solver (default = 1/2 Å)
[DEBUG  ] |  nbuffer              = 0.0                                            # Sets how far away (in grid units) the boundary of the finite difference grid is away from the solute surface
[DEBUG  ] |  nfocus               = 2                                              # Electrostatic focusing calculation
[DEBUG  ] |  fscale               = 8                                              # Set the ratio between the coarse and fine grid spacings in an electrostatic focussing calculation
[DEBUG  ] |  npbgrid              = 1                                              # Sets how often the finite-difference grid is regenerated
[DEBUG  ] |  bcopt                = 5                                              # Boundary condition option
[DEBUG  ] |  eneopt               = 2                                              # Compute electrostatic energy and forces
[DEBUG  ] |  frcopt               = 0                                              # Output for computing electrostatic forces
[DEBUG  ] |  scalec               = 0                                              # Option to compute reaction field energy and forces
[DEBUG  ] |  cutfd                = 5.0                                            # Cutoff for finite-difference interactions
[DEBUG  ] |  cutnb                = 0.0                                            # Cutoff for nonbonded interations
[DEBUG  ] |  nsnba                = 1                                              # Sets how often atom-based pairlist is generated
[DEBUG  ] |  decompopt            = 2                                              # Option to select different decomposition schemes when INP = 2
[DEBUG  ] |  use_rmin             = 1                                              # The option to set up van der Waals radii
[DEBUG  ] |  sprob                = 0.557                                          # Solvent probe radius for SASA used to compute the dispersion term
[DEBUG  ] |  vprob                = 1.3                                            # Solvent probe radius for molecular volume (the volume enclosed by SASA)
[DEBUG  ] |  rhow_effect          = 1.129                                          # Effective water density used in the non-polar dispersion term calculation
[DEBUG  ] |  use_sav              = 1                                              # Use molecular volume (the volume enclosed by SASA) for cavity term calculation
[DEBUG  ] |  cavity_surften       = 0.0378                                         # Surface tension
[DEBUG  ] |  cavity_offset        = -0.5692                                        # Offset for nonpolar solvation calc
[DEBUG  ] |  maxsph               = 400                                            # Approximate number of dots to represent the maximum atomic solvent accessible surface
[DEBUG  ] |  maxarcdot            = 1500                                           # Number of dots used to store arc dots per atom
[DEBUG  ] |  npbverb              = 0                                              # Option to turn on verbose mode
[DEBUG  ] |/
[DEBUG  ] |
[DEBUG  ] |# Decomposition namelist variables
[DEBUG  ] |&decomposition
[DEBUG  ] |  idecomp              = 1                                              # Which type of decomposition analysis to do
[DEBUG  ] |  dec_verbose          = 0                                              # Control energy terms are printed to the output
[DEBUG  ] |  print_res            = "within 6"                                     # Which residues to print decomposition data for
[DEBUG  ] |  csv_format           = 1                                              # Write decomposition data in CSV format
[DEBUG  ] |/
[DEBUG  ] |
[DEBUG  ] |--------------------------------------------------------------
[DEBUG  ] 
[INFO   ] Checking mmpbsa.in input file...
[INFO   ] Checking mmpbsa.in input file...Done.

[INFO   ] Checking external programs...
[INFO   ] cpptraj found! Using /home/romi/anaconda3/envs/gmxMMPBSA/bin/cpptraj
[INFO   ] tleap found! Using /home/romi/anaconda3/envs/gmxMMPBSA/bin/tleap
[INFO   ] parmchk2 found! Using /home/romi/anaconda3/envs/gmxMMPBSA/bin/parmchk2
[INFO   ] sander found! Using /home/romi/anaconda3/envs/gmxMMPBSA/bin/sander
[INFO   ] Using GROMACS version > 5.x.x!
[INFO   ] gmx found! Using /usr/local/gromacs/bin/gmx
[INFO   ] Checking external programs...Done.

[INFO   ] Building AMBER topologies from GROMACS files...
[INFO   ] Get PDB files from GROMACS structures files...
[INFO   ] Making gmx_MMPBSA index for complex...
[DEBUG  ] Running command: echo name 18 GMXMMPBSA_REC\n name 19 GMXMMPBSA_LIG\n  18 | 19\n q\n | /usr/local/gromacs/bin/gmx make_ndx -n index.ndx -o _GMXMMPBSA_COM_index.ndx
[DEBUG  ]                      :-) GROMACS - gmx make_ndx, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx make_ndx, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Passport/November 3rd 2023/6umt-rs1-310k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx make_ndx -n index.ndx -o _GMXMMPBSA_COM_index.ndx
[DEBUG  ] 
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "Those who cannot remember the past are condemned to compute it." (Steve Pinker)
[DEBUG  ] 
[DEBUG  ] Going to read 1 old index file(s)
[DEBUG  ] Counted atom numbers up to 65298 in index file
[DEBUG  ] 
[DEBUG  ]   0 System              : 65298 atoms
[DEBUG  ]   1 Protein             :  4113 atoms
[DEBUG  ]   2 Protein-H           :  2065 atoms
[DEBUG  ]   3 C-alpha             :   263 atoms
[DEBUG  ]   4 Backbone            :   789 atoms
[DEBUG  ]   5 MainChain           :  1050 atoms
[DEBUG  ]   6 MainChain+Cb        :  1299 atoms
[DEBUG  ]   7 MainChain+H         :  1303 atoms
[DEBUG  ]   8 SideChain           :  2810 atoms
[DEBUG  ]   9 SideChain-H         :  1015 atoms
[DEBUG  ]  10 Prot-Masses         :  4113 atoms
[DEBUG  ]  11 non-Protein         : 61185 atoms
[DEBUG  ]  12 Water               : 61179 atoms
[DEBUG  ]  13 SOL                 : 61179 atoms
[DEBUG  ]  14 non-Water           :  4119 atoms
[DEBUG  ]  15 Ion                 :     6 atoms
[DEBUG  ]  16 NA                  :     6 atoms
[DEBUG  ]  17 Water_and_ions      : 61185 atoms
[DEBUG  ]  18 Chain_A             :  2142 atoms
[DEBUG  ]  19 Chain_B             :  1971 atoms
[DEBUG  ] 
[DEBUG  ]  nr : group      '!': not  'name' nr name   'splitch' nr    Enter: list groups
[DEBUG  ]  'a': atom       '&': and  'del' nr         'splitres' nr   'l': list residues
[DEBUG  ]  't': atom type  '|': or   'keep' nr        'splitat' nr    'h': help
[DEBUG  ]  'r': residue              'res' nr         'chain' char
[DEBUG  ]  "name": group             'case': case sensitive           'q': save and quit
[DEBUG  ]  'ri': residue index
[DEBUG  ] 
[DEBUG  ] > 
[DEBUG  ] 
[DEBUG  ] > 
[DEBUG  ] 
[DEBUG  ] > 
[DEBUG  ] Copied index group 18 'GMXMMPBSA_REC'
[DEBUG  ] Copied index group 19 'GMXMMPBSA_LIG'
[DEBUG  ] Merged two groups with OR: 2142 1971 -> 4113
[DEBUG  ] 
[DEBUG  ]  20 GMXMMPBSA_REC_GMXMMPBSA_LIG:  4113 atoms
[DEBUG  ] 
[DEBUG  ] > 
[INFO   ] Normal Complex: Saving group 18_19 in _GMXMMPBSA_COM_index.ndx file as _GMXMMPBSA_COM.pdb
[DEBUG  ] Running command: echo GMXMMPBSA_REC_GMXMMPBSA_LIG | /usr/local/gromacs/bin/gmx trjconv -f system_reduced.xtc -s md_0_10.tpr -o _GMXMMPBSA_COM.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ]                      :-) GROMACS - gmx trjconv, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx trjconv, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Passport/November 3rd 2023/6umt-rs1-310k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx trjconv -f system_reduced.xtc -s md_0_10.tpr -o _GMXMMPBSA_COM.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ] 
[DEBUG  ] Will write pdb: Protein data bank file
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Group     0 (         System) has 65298 elements
[DEBUG  ] Group     1 (        Protein) has  4113 elements
[DEBUG  ] Group     2 (      Protein-H) has  2065 elements
[DEBUG  ] Group     3 (        C-alpha) has   263 elements
[DEBUG  ] Group     4 (       Backbone) has   789 elements
[DEBUG  ] Group     5 (      MainChain) has  1050 elements
[DEBUG  ] Group     6 (   MainChain+Cb) has  1299 elements
[DEBUG  ] Group     7 (    MainChain+H) has  1303 elements
[DEBUG  ] Group     8 (      SideChain) has  2810 elements
[DEBUG  ] Group     9 (    SideChain-H) has  1015 elements
[DEBUG  ] Group    10 (    Prot-Masses) has  4113 elements
[DEBUG  ] Group    11 (    non-Protein) has 61185 elements
[DEBUG  ] Group    12 (          Water) has 61179 elements
[DEBUG  ] Group    13 (            SOL) has 61179 elements
[DEBUG  ] Group    14 (      non-Water) has  4119 elements
[DEBUG  ] Group    15 (            Ion) has     6 elements
[DEBUG  ] Group    16 (             NA) has     6 elements
[DEBUG  ] Group    17 ( Water_and_ions) has 61185 elements
[DEBUG  ] Group    18 (  GMXMMPBSA_REC) has  2142 elements
[DEBUG  ] Group    19 (  GMXMMPBSA_LIG) has  1971 elements
[DEBUG  ] Group    20 (GMXMMPBSA_REC_GMXMMPBSA_LIG) has  4113 elements
[DEBUG  ] Select a group: Reading frame       0 time    0.000   
[DEBUG  ] Precision of system_reduced.xtc is 0.001 (nm)
[DEBUG  ] Reading frame       1 time  100.000   Reading frame       0 time    0.000   
[DEBUG  ] Dumping frame at t= 0 ps
[DEBUG  ]  ->  frame      0 time    0.000      Reading frame       1 time  100.000   
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "Those who cannot remember the past are condemned to compute it." (Steve Pinker)
[DEBUG  ] 
[DEBUG  ] Note that major changes are planned in future for trjconv, to improve usability and utility.
[DEBUG  ] Select group for output
[DEBUG  ] Selected 20: 'GMXMMPBSA_REC_GMXMMPBSA_LIG'
[INFO   ] No receptor structure file was defined. Using ST approach...
[INFO   ] Using receptor structure from complex to generate AMBER topology
[INFO   ] Normal Receptor: Saving group 18 in _GMXMMPBSA_COM_index.ndx file as _GMXMMPBSA_REC.pdb
[DEBUG  ] Running command: echo 18 | /usr/local/gromacs/bin/gmx trjconv -f system_reduced.xtc -s md_0_10.tpr -o _GMXMMPBSA_REC.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ]                      :-) GROMACS - gmx trjconv, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx trjconv, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Passport/November 3rd 2023/6umt-rs1-310k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx trjconv -f system_reduced.xtc -s md_0_10.tpr -o _GMXMMPBSA_REC.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ] 
[DEBUG  ] Will write pdb: Protein data bank file
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Group     0 (         System) has 65298 elements
[DEBUG  ] Group     1 (        Protein) has  4113 elements
[DEBUG  ] Group     2 (      Protein-H) has  2065 elements
[DEBUG  ] Group     3 (        C-alpha) has   263 elements
[DEBUG  ] Group     4 (       Backbone) has   789 elements
[DEBUG  ] Group     5 (      MainChain) has  1050 elements
[DEBUG  ] Group     6 (   MainChain+Cb) has  1299 elements
[DEBUG  ] Group     7 (    MainChain+H) has  1303 elements
[DEBUG  ] Group     8 (      SideChain) has  2810 elements
[DEBUG  ] Group     9 (    SideChain-H) has  1015 elements
[DEBUG  ] Group    10 (    Prot-Masses) has  4113 elements
[DEBUG  ] Group    11 (    non-Protein) has 61185 elements
[DEBUG  ] Group    12 (          Water) has 61179 elements
[DEBUG  ] Group    13 (            SOL) has 61179 elements
[DEBUG  ] Group    14 (      non-Water) has  4119 elements
[DEBUG  ] Group    15 (            Ion) has     6 elements
[DEBUG  ] Group    16 (             NA) has     6 elements
[DEBUG  ] Group    17 ( Water_and_ions) has 61185 elements
[DEBUG  ] Group    18 (  GMXMMPBSA_REC) has  2142 elements
[DEBUG  ] Group    19 (  GMXMMPBSA_LIG) has  1971 elements
[DEBUG  ] Group    20 (GMXMMPBSA_REC_GMXMMPBSA_LIG) has  4113 elements
[DEBUG  ] Select a group: Reading frame       0 time    0.000   
[DEBUG  ] Precision of system_reduced.xtc is 0.001 (nm)
[DEBUG  ] Reading frame       1 time  100.000   Reading frame       0 time    0.000   
[DEBUG  ] Dumping frame at t= 0 ps
[DEBUG  ]  ->  frame      0 time    0.000      Reading frame       1 time  100.000   
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "Those who cannot remember the past are condemned to compute it." (Steve Pinker)
[DEBUG  ] 
[DEBUG  ] Note that major changes are planned in future for trjconv, to improve usability and utility.
[DEBUG  ] Select group for output
[DEBUG  ] Selected 18: 'GMXMMPBSA_REC'
[INFO   ] No ligand structure file was defined. Using ST approach...
[INFO   ] Using ligand structure from complex to generate AMBER topology
[INFO   ] Normal Ligand: Saving group 19 in _GMXMMPBSA_COM_index.ndx file as _GMXMMPBSA_LIG.pdb
[DEBUG  ] Running command: echo 19 | /usr/local/gromacs/bin/gmx trjconv -f system_reduced.xtc -s md_0_10.tpr -o _GMXMMPBSA_LIG.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ]                      :-) GROMACS - gmx trjconv, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx trjconv, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Passport/November 3rd 2023/6umt-rs1-310k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx trjconv -f system_reduced.xtc -s md_0_10.tpr -o _GMXMMPBSA_LIG.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ] 
[DEBUG  ] Will write pdb: Protein data bank file
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Group     0 (         System) has 65298 elements
[DEBUG  ] Group     1 (        Protein) has  4113 elements
[DEBUG  ] Group     2 (      Protein-H) has  2065 elements
[DEBUG  ] Group     3 (        C-alpha) has   263 elements
[DEBUG  ] Group     4 (       Backbone) has   789 elements
[DEBUG  ] Group     5 (      MainChain) has  1050 elements
[DEBUG  ] Group     6 (   MainChain+Cb) has  1299 elements
[DEBUG  ] Group     7 (    MainChain+H) has  1303 elements
[DEBUG  ] Group     8 (      SideChain) has  2810 elements
[DEBUG  ] Group     9 (    SideChain-H) has  1015 elements
[DEBUG  ] Group    10 (    Prot-Masses) has  4113 elements
[DEBUG  ] Group    11 (    non-Protein) has 61185 elements
[DEBUG  ] Group    12 (          Water) has 61179 elements
[DEBUG  ] Group    13 (            SOL) has 61179 elements
[DEBUG  ] Group    14 (      non-Water) has  4119 elements
[DEBUG  ] Group    15 (            Ion) has     6 elements
[DEBUG  ] Group    16 (             NA) has     6 elements
[DEBUG  ] Group    17 ( Water_and_ions) has 61185 elements
[DEBUG  ] Group    18 (  GMXMMPBSA_REC) has  2142 elements
[DEBUG  ] Group    19 (  GMXMMPBSA_LIG) has  1971 elements
[DEBUG  ] Group    20 (GMXMMPBSA_REC_GMXMMPBSA_LIG) has  4113 elements
[DEBUG  ] Select a group: Reading frame       0 time    0.000   
[DEBUG  ] Precision of system_reduced.xtc is 0.001 (nm)
[DEBUG  ] Reading frame       1 time  100.000   Reading frame       0 time    0.000   
[DEBUG  ] Dumping frame at t= 0 ps
[DEBUG  ]  ->  frame      0 time    0.000      Reading frame       1 time  100.000   
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "Those who cannot remember the past are condemned to compute it." (Steve Pinker)
[DEBUG  ] 
[DEBUG  ] Note that major changes are planned in future for trjconv, to improve usability and utility.
[DEBUG  ] Select group for output
[DEBUG  ] Selected 19: 'GMXMMPBSA_LIG'
[INFO   ] Checking the structures consistency...
[INFO   ] 
[INFO   ] Using topology conversion. Setting radiopt = 0...
[INFO   ] Building Normal Complex Amber topology...
[INFO   ] Detected CHARMM force field topology format...
[INFO   ] Assigning PBRadii mbondi2 to Complex...
[INFO   ] Writing Normal Complex AMBER topology...
[INFO   ] No Receptor topology file was defined. Using ST approach...
[INFO   ] Building AMBER Receptor topology from Complex...
[INFO   ] Assigning PBRadii mbondi2 to Receptor...
[INFO   ] Writing Normal Receptor AMBER topology...
[INFO   ] No Ligand topology file was defined. Using ST approach...
[INFO   ] Building AMBER Ligand topology from Complex...
[INFO   ] Assigning PBRadii mbondi2 to Ligand...
[INFO   ] Writing Normal Ligand AMBER topology...
[INFO   ] Selecting residues by distance (6 Å) between receptor and ligand for decomposition analysis...
[INFO   ] Selected 90 residues:
R:A:GLU:48  R:A:SER:49  R:A:PHE:50  R:A:VAL:51  R:A:LEU:52  R:A:ASN:53  R:A:TYR:55  R:A:MET:57  R:A:SER:58  R:A:SER:60
R:A:GLY:61  R:A:GLN:62  R:A:PRO:63  R:A:ASP:64  R:A:LYS:65  R:A:LEU:66  R:A:ALA:67  R:A:ALA:68  R:A:PRO:70  R:A:GLU:71
R:A:ASP:72  R:A:ARG:73  R:A:SER:74  R:A:GLN:75  R:A:PRO:76  R:A:GLY:77  R:A:GLN:78  R:A:ASP:79  R:A:SER:80  R:A:ARG:81
R:A:PHE:82  R:A:LEU:109 R:A:CYS:110 R:A:GLY:111 R:A:ALA:112 R:A:ILE:113 R:A:LEU:115 R:A:PRO:117 R:A:LYS:118 R:A:VAL:119
R:A:GLN:120 R:A:ILE:121 R:A:LYS:122 R:A:GLU:123 R:A:ARG:126 L:B:LEU:152 L:B:GLN:153 L:B:LEU:154 L:B:HIS:155 L:B:GLN:156
L:B:ILE:157 L:B:ALA:158 L:B:ALA:159 L:B:LEU:160 L:B:PHE:161 L:B:THR:162 L:B:VAL:163 L:B:THR:164 L:B:VAL:165 L:B:GLU:168
L:B:ASP:185 L:B:THR:186 L:B:GLY:187 L:B:SER:188 L:B:THR:196 L:B:ALA:197 L:B:SER:198 L:B:GLN:200 L:B:ASP:205 L:B:THR:206
L:B:SER:207 L:B:PRO:208 L:B:HIS:209 L:B:ARG:210 L:B:GLU:211 L:B:GLN:239 L:B:GLN:241 L:B:CYS:242 L:B:ILE:243 L:B:ILE:244
L:B:ILE:245 L:B:TYR:246 L:B:GLY:247 L:B:VAL:248 L:B:ALA:249 L:B:TRP:250 L:B:ASP:251 L:B:TYR:252 L:B:LYS:253 L:B:TYR:254

[INFO   ] Cleaning normal complex trajectories...
[DEBUG  ] Running command: echo GMXMMPBSA_REC_GMXMMPBSA_LIG | /usr/local/gromacs/bin/gmx trjconv -f system_reduced.xtc -s md_0_10.tpr -o COM_traj_0.xtc -n _GMXMMPBSA_COM_index.ndx
[DEBUG  ]                      :-) GROMACS - gmx trjconv, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx trjconv, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Passport/November 3rd 2023/6umt-rs1-310k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx trjconv -f system_reduced.xtc -s md_0_10.tpr -o COM_traj_0.xtc -n _GMXMMPBSA_COM_index.ndx
[DEBUG  ] 
[DEBUG  ] Will write xtc: Compressed trajectory (portable xdr format): xtc
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Group     0 (         System) has 65298 elements
[DEBUG  ] Group     1 (        Protein) has  4113 elements
[DEBUG  ] Group     2 (      Protein-H) has  2065 elements
[DEBUG  ] Group     3 (        C-alpha) has   263 elements
[DEBUG  ] Group     4 (       Backbone) has   789 elements
[DEBUG  ] Group     5 (      MainChain) has  1050 elements
[DEBUG  ] Group     6 (   MainChain+Cb) has  1299 elements
[DEBUG  ] Group     7 (    MainChain+H) has  1303 elements
[DEBUG  ] Group     8 (      SideChain) has  2810 elements
[DEBUG  ] Group     9 (    SideChain-H) has  1015 elements
[DEBUG  ] Group    10 (    Prot-Masses) has  4113 elements
[DEBUG  ] Group    11 (    non-Protein) has 61185 elements
[DEBUG  ] Group    12 (          Water) has 61179 elements
[DEBUG  ] Group    13 (            SOL) has 61179 elements
[DEBUG  ] Group    14 (      non-Water) has  4119 elements
[DEBUG  ] Group    15 (            Ion) has     6 elements
[DEBUG  ] Group    16 (             NA) has     6 elements
[DEBUG  ] Group    17 ( Water_and_ions) has 61185 elements
[DEBUG  ] Group    18 (  GMXMMPBSA_REC) has  2142 elements
[DEBUG  ] Group    19 (  GMXMMPBSA_LIG) has  1971 elements
[DEBUG  ] Group    20 (GMXMMPBSA_REC_GMXMMPBSA_LIG) has  4113 elements
[DEBUG  ] Select a group: Reading frame       0 time    0.000   
[DEBUG  ] Precision of system_reduced.xtc is 0.001 (nm)
[DEBUG  ] Using output precision of 0.001 (nm)
[DEBUG  ]  ->  frame      0 time    0.000      Reading frame       1 time  100.000    ->  frame      1 time  100.000      Reading frame       2 time  200.000    ->  frame      2 time  200.000      Reading frame       3 time  300.000    ->  frame      3 time  300.000      Reading frame       4 time  400.000    ->  frame      4 time  400.000      Reading frame       5 time  500.000    ->  frame      5 time  500.000      Reading frame       6 time  600.000    ->  frame      6 time  600.000      Reading frame       7 time  700.000    ->  frame      7 time  700.000      Reading frame       8 time  800.000    ->  frame      8 time  800.000      Reading frame       9 time  900.000    ->  frame      9 time  900.000      Reading frame      10 time 1000.000    ->  frame     10 time 1000.000      Reading frame      11 time 1100.000   Reading frame      12 time 1200.000   Reading frame      13 time 1300.000   Reading frame      14 time 1400.000   Reading frame      15 time 1500.000   Reading frame      16 time 1600.000   Reading frame      17 time 1700.000   Reading frame      18 time 1800.000   Reading frame      19 time 1900.000   Reading frame      20 time 2000.000    ->  frame     20 time 2000.000      Reading frame      30 time 3000.000    ->  frame     30 time 3000.000      Reading frame      40 time 4000.000    ->  frame     40 time 4000.000      Reading frame      50 time 5000.000    ->  frame     50 time 5000.000      Reading frame      60 time 6000.000    ->  frame     60 time 6000.000      Reading frame      70 time 7000.000    ->  frame     70 time 7000.000      Reading frame      80 time 8000.000    ->  frame     80 time 8000.000      Reading frame      90 time 9000.000    ->  frame     90 time 9000.000      Reading frame     100 time 10000.000    ->  frame    100 time 10000.000      Reading frame     110 time 11000.000    ->  frame    110 time 11000.000      Reading frame     120 time 12000.000    ->  frame    120 time 12000.000      Reading frame     130 time 13000.000    ->  frame    130 time 13000.000      Reading frame     140 time 14000.000    ->  frame    140 time 14000.000      Reading frame     150 time 15000.000    ->  frame    150 time 15000.000      Reading frame     160 time 16000.000    ->  frame    160 time 16000.000      Reading frame     170 time 17000.000    ->  frame    170 time 17000.000      Reading frame     180 time 18000.000    ->  frame    180 time 18000.000      Reading frame     190 time 19000.000    ->  frame    190 time 19000.000      Reading frame     200 time 20000.000    ->  frame    200 time 20000.000       ->  frame    210 time 21000.000       ->  frame    220 time 22000.000       ->  frame    230 time 23000.000       ->  frame    240 time 24000.000       ->  frame    250 time 25000.000       ->  frame    260 time 26000.000       ->  frame    270 time 27000.000       ->  frame    280 time 28000.000       ->  frame    290 time 29000.000      Reading frame     300 time 30000.000    ->  frame    300 time 30000.000       ->  frame    310 time 31000.000       ->  frame    320 time 32000.000       ->  frame    330 time 33000.000       ->  frame    340 time 34000.000       ->  frame    350 time 35000.000       ->  frame    360 time 36000.000       ->  frame    370 time 37000.000       ->  frame    380 time 38000.000       ->  frame    390 time 39000.000      Reading frame     400 time 40000.000    ->  frame    400 time 40000.000       ->  frame    410 time 41000.000       ->  frame    420 time 42000.000       ->  frame    430 time 43000.000       ->  frame    440 time 44000.000       ->  frame    450 time 45000.000       ->  frame    460 time 46000.000       ->  frame    470 time 47000.000       ->  frame    480 time 48000.000       ->  frame    490 time 49000.000      Reading frame     500 time 50000.000    ->  frame    500 time 50000.000       ->  frame    510 time 51000.000       ->  frame    520 time 52000.000       ->  frame    530 time 53000.000       ->  frame    540 time 54000.000       ->  frame    550 time 55000.000       ->  frame    560 time 56000.000       ->  frame    570 time 57000.000       ->  frame    580 time 58000.000       ->  frame    590 time 59000.000      Reading frame     600 time 60000.000    ->  frame    600 time 60000.000       ->  frame    610 time 61000.000       ->  frame    620 time 62000.000       ->  frame    630 time 63000.000       ->  frame    640 time 64000.000       ->  frame    650 time 65000.000       ->  frame    660 time 66000.000       ->  frame    670 time 67000.000       ->  frame    680 time 68000.000       ->  frame    690 time 69000.000      Reading frame     700 time 70000.000    ->  frame    700 time 70000.000       ->  frame    710 time 71000.000       ->  frame    720 time 72000.000       ->  frame    730 time 73000.000       ->  frame    740 time 74000.000       ->  frame    750 time 75000.000       ->  frame    760 time 76000.000       ->  frame    770 time 77000.000       ->  frame    780 time 78000.000       ->  frame    790 time 79000.000      Reading frame     800 time 80000.000    ->  frame    800 time 80000.000       ->  frame    810 time 81000.000       ->  frame    820 time 82000.000       ->  frame    830 time 83000.000       ->  frame    840 time 84000.000       ->  frame    850 time 85000.000       ->  frame    860 time 86000.000       ->  frame    870 time 87000.000       ->  frame    880 time 88000.000       ->  frame    890 time 89000.000      Reading frame     900 time 90000.000    ->  frame    900 time 90000.000       ->  frame    910 time 91000.000       ->  frame    920 time 92000.000       ->  frame    930 time 93000.000       ->  frame    940 time 94000.000       ->  frame    950 time 95000.000       ->  frame    960 time 96000.000       ->  frame    970 time 97000.000       ->  frame    980 time 98000.000       ->  frame    990 time 99000.000      Reading frame    1000 time 100000.000    ->  frame   1000 time 100000.000       ->  frame   1010 time 101000.000       ->  frame   1020 time 102000.000       ->  frame   1030 time 103000.000       ->  frame   1040 time 104000.000       ->  frame   1050 time 105000.000       ->  frame   1060 time 106000.000       ->  frame   1070 time 107000.000       ->  frame   1080 time 108000.000       ->  frame   1090 time 109000.000      Reading frame    1100 time 110000.000    ->  frame   1100 time 110000.000       ->  frame   1110 time 111000.000       ->  frame   1120 time 112000.000       ->  frame   1130 time 113000.000       ->  frame   1140 time 114000.000       ->  frame   1150 time 115000.000       ->  frame   1160 time 116000.000       ->  frame   1170 time 117000.000       ->  frame   1180 time 118000.000       ->  frame   1190 time 119000.000      Reading frame    1200 time 120000.000    ->  frame   1200 time 120000.000       ->  frame   1210 time 121000.000       ->  frame   1220 time 122000.000       ->  frame   1230 time 123000.000       ->  frame   1240 time 124000.000       ->  frame   1250 time 125000.000       ->  frame   1260 time 126000.000       ->  frame   1270 time 127000.000       ->  frame   1280 time 128000.000       ->  frame   1290 time 129000.000      Reading frame    1300 time 130000.000    ->  frame   1300 time 130000.000       ->  frame   1310 time 131000.000       ->  frame   1320 time 132000.000       ->  frame   1330 time 133000.000       ->  frame   1340 time 134000.000       ->  frame   1350 time 135000.000       ->  frame   1360 time 136000.000       ->  frame   1370 time 137000.000       ->  frame   1380 time 138000.000       ->  frame   1390 time 139000.000      Reading frame    1400 time 140000.000    ->  frame   1400 time 140000.000       ->  frame   1410 time 141000.000       ->  frame   1420 time 142000.000       ->  frame   1430 time 143000.000       ->  frame   1440 time 144000.000       ->  frame   1450 time 145000.000       ->  frame   1460 time 146000.000       ->  frame   1470 time 147000.000       ->  frame   1480 time 148000.000       ->  frame   1490 time 149000.000      Reading frame    1500 time 150000.000    ->  frame   1500 time 150000.000       ->  frame   1510 time 151000.000       ->  frame   1520 time 152000.000       ->  frame   1530 time 153000.000       ->  frame   1540 time 154000.000       ->  frame   1550 time 155000.000       ->  frame   1560 time 156000.000       ->  frame   1570 time 157000.000       ->  frame   1580 time 158000.000       ->  frame   1590 time 159000.000      Reading frame    1600 time 160000.000    ->  frame   1600 time 160000.000       ->  frame   1610 time 161000.000       ->  frame   1620 time 162000.000       ->  frame   1630 time 163000.000       ->  frame   1640 time 164000.000       ->  frame   1650 time 165000.000       ->  frame   1660 time 166000.000       ->  frame   1670 time 167000.000       ->  frame   1680 time 168000.000       ->  frame   1690 time 169000.000      Reading frame    1700 time 170000.000    ->  frame   1700 time 170000.000       ->  frame   1710 time 171000.000       ->  frame   1720 time 172000.000       ->  frame   1730 time 173000.000       ->  frame   1740 time 174000.000       ->  frame   1750 time 175000.000       ->  frame   1760 time 176000.000       ->  frame   1770 time 177000.000       ->  frame   1780 time 178000.000       ->  frame   1790 time 179000.000      Reading frame    1800 time 180000.000    ->  frame   1800 time 180000.000       ->  frame   1810 time 181000.000       ->  frame   1820 time 182000.000       ->  frame   1830 time 183000.000       ->  frame   1840 time 184000.000       ->  frame   1850 time 185000.000       ->  frame   1860 time 186000.000       ->  frame   1870 time 187000.000       ->  frame   1880 time 188000.000       ->  frame   1890 time 189000.000      Reading frame    1900 time 190000.000    ->  frame   1900 time 190000.000       ->  frame   1910 time 191000.000       ->  frame   1920 time 192000.000       ->  frame   1930 time 193000.000       ->  frame   1940 time 194000.000       ->  frame   1950 time 195000.000       ->  frame   1960 time 196000.000       ->  frame   1970 time 197000.000       ->  frame   1980 time 198000.000       ->  frame   1990 time 199000.000      Reading frame    2000 time 200000.000    ->  frame   2000 time 200000.000       ->  frame   2010 time 201000.000       ->  frame   2020 time 202000.000       ->  frame   2030 time 203000.000       ->  frame   2040 time 204000.000       ->  frame   2050 time 205000.000       ->  frame   2060 time 206000.000       ->  frame   2070 time 207000.000       ->  frame   2080 time 208000.000       ->  frame   2090 time 209000.000       ->  frame   2100 time 210000.000       ->  frame   2110 time 211000.000       ->  frame   2120 time 212000.000       ->  frame   2130 time 213000.000       ->  frame   2140 time 214000.000       ->  frame   2150 time 215000.000       ->  frame   2160 time 216000.000       ->  frame   2170 time 217000.000       ->  frame   2180 time 218000.000       ->  frame   2190 time 219000.000       ->  frame   2200 time 220000.000       ->  frame   2210 time 221000.000       ->  frame   2220 time 222000.000       ->  frame   2230 time 223000.000       ->  frame   2240 time 224000.000       ->  frame   2250 time 225000.000       ->  frame   2260 time 226000.000       ->  frame   2270 time 227000.000       ->  frame   2280 time 228000.000       ->  frame   2290 time 229000.000       ->  frame   2300 time 230000.000       ->  frame   2310 time 231000.000       ->  frame   2320 time 232000.000       ->  frame   2330 time 233000.000       ->  frame   2340 time 234000.000       ->  frame   2350 time 235000.000       ->  frame   2360 time 236000.000       ->  frame   2370 time 237000.000       ->  frame   2380 time 238000.000       ->  frame   2390 time 239000.000       ->  frame   2400 time 240000.000       ->  frame   2410 time 241000.000       ->  frame   2420 time 242000.000       ->  frame   2430 time 243000.000       ->  frame   2440 time 244000.000       ->  frame   2450 time 245000.000       ->  frame   2460 time 246000.000       ->  frame   2470 time 247000.000       ->  frame   2480 time 248000.000       ->  frame   2490 time 249000.000       ->  frame   2500 time 250000.000       ->  frame   2510 time 251000.000       ->  frame   2520 time 252000.000       ->  frame   2530 time 253000.000       ->  frame   2540 time 254000.000       ->  frame   2550 time 255000.000       ->  frame   2560 time 256000.000       ->  frame   2570 time 257000.000       ->  frame   2580 time 258000.000       ->  frame   2590 time 259000.000       ->  frame   2600 time 260000.000       ->  frame   2610 time 261000.000       ->  frame   2620 time 262000.000       ->  frame   2630 time 263000.000       ->  frame   2640 time 264000.000       ->  frame   2650 time 265000.000       ->  frame   2660 time 266000.000       ->  frame   2670 time 267000.000       ->  frame   2680 time 268000.000       ->  frame   2690 time 269000.000       ->  frame   2700 time 270000.000       ->  frame   2710 time 271000.000       ->  frame   2720 time 272000.000       ->  frame   2730 time 273000.000       ->  frame   2740 time 274000.000       ->  frame   2750 time 275000.000       ->  frame   2760 time 276000.000       ->  frame   2770 time 277000.000       ->  frame   2780 time 278000.000       ->  frame   2790 time 279000.000       ->  frame   2800 time 280000.000       ->  frame   2810 time 281000.000       ->  frame   2820 time 282000.000       ->  frame   2830 time 283000.000       ->  frame   2840 time 284000.000       ->  frame   2850 time 285000.000       ->  frame   2860 time 286000.000       ->  frame   2870 time 287000.000       ->  frame   2880 time 288000.000       ->  frame   2890 time 289000.000       ->  frame   2900 time 290000.000       ->  frame   2910 time 291000.000       ->  frame   2920 time 292000.000       ->  frame   2930 time 293000.000       ->  frame   2940 time 294000.000       ->  frame   2950 time 295000.000       ->  frame   2960 time 296000.000       ->  frame   2970 time 297000.000       ->  frame   2980 time 298000.000       ->  frame   2990 time 299000.000      Reading frame    3000 time 300000.000    ->  frame   3000 time 300000.000       ->  frame   3010 time 301000.000       ->  frame   3020 time 302000.000       ->  frame   3030 time 303000.000       ->  frame   3040 time 304000.000       ->  frame   3050 time 305000.000       ->  frame   3060 time 306000.000       ->  frame   3070 time 307000.000       ->  frame   3080 time 308000.000       ->  frame   3090 time 309000.000       ->  frame   3100 time 310000.000       ->  frame   3110 time 311000.000       ->  frame   3120 time 312000.000       ->  frame   3130 time 313000.000       ->  frame   3140 time 314000.000       ->  frame   3150 time 315000.000       ->  frame   3160 time 316000.000       ->  frame   3170 time 317000.000       ->  frame   3180 time 318000.000       ->  frame   3190 time 319000.000       ->  frame   3200 time 320000.000       ->  frame   3210 time 321000.000       ->  frame   3220 time 322000.000       ->  frame   3230 time 323000.000       ->  frame   3240 time 324000.000       ->  frame   3250 time 325000.000       ->  frame   3260 time 326000.000       ->  frame   3270 time 327000.000       ->  frame   3280 time 328000.000       ->  frame   3290 time 329000.000       ->  frame   3300 time 330000.000       ->  frame   3310 time 331000.000       ->  frame   3320 time 332000.000       ->  frame   3330 time 333000.000       ->  frame   3340 time 334000.000       ->  frame   3350 time 335000.000       ->  frame   3360 time 336000.000       ->  frame   3370 time 337000.000       ->  frame   3380 time 338000.000       ->  frame   3390 time 339000.000       ->  frame   3400 time 340000.000       ->  frame   3410 time 341000.000       ->  frame   3420 time 342000.000       ->  frame   3430 time 343000.000       ->  frame   3440 time 344000.000       ->  frame   3450 time 345000.000       ->  frame   3460 time 346000.000       ->  frame   3470 time 347000.000       ->  frame   3480 time 348000.000       ->  frame   3490 time 349000.000       ->  frame   3500 time 350000.000       ->  frame   3510 time 351000.000       ->  frame   3520 time 352000.000       ->  frame   3530 time 353000.000       ->  frame   3540 time 354000.000       ->  frame   3550 time 355000.000       ->  frame   3560 time 356000.000       ->  frame   3570 time 357000.000       ->  frame   3580 time 358000.000       ->  frame   3590 time 359000.000       ->  frame   3600 time 360000.000       ->  frame   3610 time 361000.000       ->  frame   3620 time 362000.000       ->  frame   3630 time 363000.000       ->  frame   3640 time 364000.000       ->  frame   3650 time 365000.000       ->  frame   3660 time 366000.000       ->  frame   3670 time 367000.000       ->  frame   3680 time 368000.000       ->  frame   3690 time 369000.000       ->  frame   3700 time 370000.000       ->  frame   3710 time 371000.000       ->  frame   3720 time 372000.000       ->  frame   3730 time 373000.000       ->  frame   3740 time 374000.000       ->  frame   3750 time 375000.000       ->  frame   3760 time 376000.000       ->  frame   3770 time 377000.000       ->  frame   3780 time 378000.000       ->  frame   3790 time 379000.000       ->  frame   3800 time 380000.000       ->  frame   3810 time 381000.000       ->  frame   3820 time 382000.000       ->  frame   3830 time 383000.000       ->  frame   3840 time 384000.000       ->  frame   3850 time 385000.000       ->  frame   3860 time 386000.000       ->  frame   3870 time 387000.000       ->  frame   3880 time 388000.000       ->  frame   3890 time 389000.000       ->  frame   3900 time 390000.000       ->  frame   3910 time 391000.000       ->  frame   3920 time 392000.000       ->  frame   3930 time 393000.000       ->  frame   3940 time 394000.000       ->  frame   3950 time 395000.000       ->  frame   3960 time 396000.000       ->  frame   3970 time 397000.000       ->  frame   3980 time 398000.000       ->  frame   3990 time 399000.000      Reading frame    4000 time 400000.000    ->  frame   4000 time 400000.000       ->  frame   4010 time 401000.000       ->  frame   4020 time 402000.000       ->  frame   4030 time 403000.000       ->  frame   4040 time 404000.000       ->  frame   4050 time 405000.000       ->  frame   4060 time 406000.000       ->  frame   4070 time 407000.000       ->  frame   4080 time 408000.000       ->  frame   4090 time 409000.000       ->  frame   4100 time 410000.000       ->  frame   4110 time 411000.000       ->  frame   4120 time 412000.000       ->  frame   4130 time 413000.000       ->  frame   4140 time 414000.000       ->  frame   4150 time 415000.000       ->  frame   4160 time 416000.000       ->  frame   4170 time 417000.000       ->  frame   4180 time 418000.000       ->  frame   4190 time 419000.000       ->  frame   4200 time 420000.000       ->  frame   4210 time 421000.000       ->  frame   4220 time 422000.000       ->  frame   4230 time 423000.000       ->  frame   4240 time 424000.000       ->  frame   4250 time 425000.000       ->  frame   4260 time 426000.000       ->  frame   4270 time 427000.000       ->  frame   4280 time 428000.000       ->  frame   4290 time 429000.000       ->  frame   4300 time 430000.000       ->  frame   4310 time 431000.000       ->  frame   4320 time 432000.000       ->  frame   4330 time 433000.000       ->  frame   4340 time 434000.000       ->  frame   4350 time 435000.000       ->  frame   4360 time 436000.000       ->  frame   4370 time 437000.000       ->  frame   4380 time 438000.000       ->  frame   4390 time 439000.000       ->  frame   4400 time 440000.000       ->  frame   4410 time 441000.000       ->  frame   4420 time 442000.000       ->  frame   4430 time 443000.000       ->  frame   4440 time 444000.000       ->  frame   4450 time 445000.000       ->  frame   4460 time 446000.000       ->  frame   4470 time 447000.000       ->  frame   4480 time 448000.000       ->  frame   4490 time 449000.000       ->  frame   4500 time 450000.000       ->  frame   4510 time 451000.000       ->  frame   4520 time 452000.000       ->  frame   4530 time 453000.000       ->  frame   4540 time 454000.000       ->  frame   4550 time 455000.000       ->  frame   4560 time 456000.000       ->  frame   4570 time 457000.000       ->  frame   4580 time 458000.000       ->  frame   4590 time 459000.000       ->  frame   4600 time 460000.000       ->  frame   4610 time 461000.000       ->  frame   4620 time 462000.000       ->  frame   4630 time 463000.000       ->  frame   4640 time 464000.000       ->  frame   4650 time 465000.000       ->  frame   4660 time 466000.000       ->  frame   4670 time 467000.000       ->  frame   4680 time 468000.000       ->  frame   4690 time 469000.000       ->  frame   4700 time 470000.000       ->  frame   4710 time 471000.000       ->  frame   4720 time 472000.000       ->  frame   4730 time 473000.000       ->  frame   4740 time 474000.000       ->  frame   4750 time 475000.000       ->  frame   4760 time 476000.000       ->  frame   4770 time 477000.000       ->  frame   4780 time 478000.000       ->  frame   4790 time 479000.000       ->  frame   4800 time 480000.000       ->  frame   4810 time 481000.000       ->  frame   4820 time 482000.000       ->  frame   4830 time 483000.000       ->  frame   4840 time 484000.000       ->  frame   4850 time 485000.000       ->  frame   4860 time 486000.000       ->  frame   4870 time 487000.000       ->  frame   4880 time 488000.000       ->  frame   4890 time 489000.000       ->  frame   4900 time 490000.000       ->  frame   4910 time 491000.000       ->  frame   4920 time 492000.000       ->  frame   4930 time 493000.000       ->  frame   4940 time 494000.000       ->  frame   4950 time 495000.000       ->  frame   4960 time 496000.000       ->  frame   4970 time 497000.000       ->  frame   4980 time 498000.000       ->  frame   4990 time 499000.000      Reading frame    5000 time 500000.000    ->  frame   5000 time 500000.000      Last frame       5000 time 500000.000   
[DEBUG  ] 
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "It is unfortunate that the authors did not make better use of all the electric power energy that went into these massive computations." (An anonymous referee)
[DEBUG  ] 
[DEBUG  ] Note that major changes are planned in future for trjconv, to improve usability and utility.
[DEBUG  ] Select group for output
[DEBUG  ] Selected 20: 'GMXMMPBSA_REC_GMXMMPBSA_LIG'
[INFO   ] Building AMBER topologies from GROMACS files... Done.

[INFO   ] Loading and checking parameter files for compatibility...
[INFO   ] Preparing trajectories for simulation...

[INFO   ] 53 frames were processed by cpptraj for use in calculation.
[INFO   ] Running calculations on normal system...
[INFO   ] Beginning GB calculations with /home/romi/anaconda3/envs/gmxMMPBSA/bin/sander
[INFO   ]   calculating complex contribution...
[INFO   ]   calculating receptor contribution...
[INFO   ]   calculating ligand contribution...
[INFO   ] Beginning PB calculations with /home/romi/anaconda3/envs/gmxMMPBSA/bin/sander
[INFO   ]   calculating complex contribution...
[INFO   ]   calculating receptor contribution...
[INFO   ]   calculating ligand contribution...
[INFO   ] Parsing results to output files...

[DEBUG  ] Creating converter from 5 to 3
[INFO   ] Timing:
[INFO   ] Total GROMACS setup time:                   0.261 min.
[INFO   ] Total AMBER setup time:                     0.009 min.
[INFO   ] Creating trajectories with cpptraj:         0.054 min.
[INFO   ] Total calculation time:                    42.752 min.
[INFO   ] Total GB calculation time:                  6.842 min.
[INFO   ] Total PB calculation time:                 35.909 min.
[INFO   ] Statistics calculation & output writing:    0.008 min.
[INFO   ] Total time taken:                          43.098 min.

[INFO   ] 
   Finalizing gmx_MMPBSA: [ERROR  ] = 0; [WARNING] = 0
   Check the gmx_MMPBSA.log file for more details...

[INFO   ] 
 Thank you for using gmx_MMPBSA. Please consider supporting gmx_MMPBSA by citing our publication:
    Valdés-Tresanco, M.S., Valdés-Tresanco, M.E., Valiente, P.A. and Moreno E. 
    gmx_MMPBSA: A New Tool to Perform End-State Free Energy Calculations with GROMACS. 
    J Chem Theory Comput., 2021, 17 (10):6281-6291. Epub 2021 Sep 29. PMID: 34586825.
    https://pubs.acs.org/doi/10.1021/acs.jctc.1c00645

Also consider citing MMPBSA.py:
    Miller III, B. R., McGee Jr., T. D., Swails, J. M. Homeyer, N. Gohlke, H. and Roitberg, A. E.
    MMPBSA.py: An Efficient Program for End-State Free Energy Calculations.
    J. Chem. Theory Comput., 2012, 8 (9) pp 3314-3321

[INFO   ] Opening gmx_MMPBSA_ana to analyze results...

[INFO   ] Finalized...
